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Abstract 

The domain of validity of standard thermodynamics and Boltzmann-Gibbs 
statistical mechanics is discussed and then formally enlarged in order to hope- 
fully cover a variety of anomalous systems. The generalization concerns 
nonextensive systems, where nonextensivity is understood in the thermody- 
namical sense. This generalization was first proposed in 1988 inspired by the 
probabilistic description of multifractal geometries, and has been intensively 
studied during this decade. In the present effort, after introducing some his- 
torical background, we briefly describe the formalism, and then exhibit the 
present status in what concerns theoretical, experimental and computational 
evidences and connections, as well as some perspectives for the future. In 
addition to these, here and there we point out various (possibly!) relevant 
questions, whose answer would certainly clarify our current understanding 
of the foundations of statistical mechanics and its thermodynamical implica- 
tions. 



I. INTRODUCTION 

A diffuse belief exists, among many physicists as well as other scientists, that Boltzmann- 
Gibbs (BG) statistical mechanics and standard thermodynamics are eternal and universal. 
It is certainly fair to say that "eternal" , in precisely the same sense that Newtonian mechan- 
ics is "eternal", they indeed are. But, again in complete analogy with Newtonian mechanics, 
we can by no means consider them as universal. Indeed, we all know that, when the involved 
velocities approach that of light, Newtonian mechanics becomes only an approximation (an 
increasingly bad one!) and reality is better described by special relativity. Analogously, 
when the involved masses are as small as say the electron mass, once again Newtonian me- 
chanics becomes but a (bad) approximation, and quantum mechanics becomes necessary to 
understand nature. Also, if the involved masses are very large, Newtonian mechanics has 
to be extended into general relativity. In these senses we certainly cannot consider New- 
tonian mechanics as being universal. I believe that the same type of considerations apply 
to standard statistical mechanics and thermodynamics. Indeed, after more than one cen- 
tury highly successful applications of the magnificent Boltzmann's connection of Clausius 
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macroscopic entropy to the theory of probabilities applied to the microscopic world , BG 
thermal statistics can (and should!) easily be considered as one of the pillars of modern 
science. However, it is unavoidable to think that, like all other products of human mind, 
this formalism must have physical restrictions, i.e., domains of applicability, out of which 
it can at best be but an approximation. It seems that BG statistics satisfactorily describe 
nature if the effective microscopic interactions are short-ianged (i.e., close spatial connec- 
tions) and the effective microscopic memory is short-remged (i.e., close time connections) and 
the boundary conditions are non(multi) fractal. Roughly speaking, the standard formalisms 
are applicable whenever (and probably only whenever) the relevant space-time (hence the 
relevant phase space) is non(multi) fractal. If this is not the case, some kind of extension 
appears to become necessary. Indeed, an everyday increasing list of physical anomalies are, 
here and there, being pointed out which defy (not to say that plainly violate!) the standard 
BG prescriptions. A nonextensive thermostatistics, which recovers the extensive, BG one 
as particular case, was proposed in 1988 which might correctly cover at least some of 
the known anomalies. Although a fair amount of what legitimately looks like being success- 
ful applications is nowadays accumulating, further verifications and deep understanding is 
needed and welcome. Computational work is highly desired since, on various grounds, the 
analytic discussion frankly appears to be untractable. Needless, of course, to say that more 
experimental and theoretical work is absolutely relevant to exhibit the applicability and ro- 
bustness of the ideas I intend to present herein. In the present contribution, I propose some 
(hopefully relevant!) questions that are right now open to such theoretical, experimental 
and computational contributions. 

Let us be more specific. As mentioned above, it is nowadays quite well known that a 
variety of physical systems exists for which the powerful (and beautiful!) BG statistical 
mechanics and standard thermodynamics present serious difficulties or anomalies, which 
can occasionally achieve the status of just plain failures. Within a long list that will be 
systematically focused on later on, we may mention at this point systems involving long- 
range interactions (e.g., d = 3 gravitation) ||, long-range microscopic memory (e.g., non- 
markovian stochastic processes, on which much remains to be known, in fact) @,[|], and, 
generally speaking, conservative (e.g., Hamiltonian) or dissipative systems which in one 
way or another involve a relevant space-time (hence, a relevant phase space) which has a 
(multi)fractal-like structure. For instance, pure-electron plasma two-dimensional turbulence 
|| , Levy anomalous diffusion |7j , granular systems |§ , phonon-electron anomalous thermal- 
ization in ion-bombarded solids ( |J and references therein), solar neutrinos [|It]], peculiar 
velocities of galaxies [11|, inverse bremsstrahlung in plasma Jl2] and black holes fl3f , to 



cite a few, clearly appear to be (in some cases), or could possibly be (in others), concrete 
examples. The present status of these and others will be discussed in Sections III, IV and 
V. 



II. FORMALISM 
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A. ENTROPY 



As an attempt to overcome at least some of these difficulties a proposal has been ad- 



vanced, one decade ago (see also |14],|15|]), which is based on a generalized entropic form, 
namely 



s q = k L ~^= lPi [5> = i; qen) , (t) 



i=i 



where k is a positive constant and W is the total number of microscopic possibilities of 
the system (for the q < case, care must be taken to exclude all those possibilities whose 
probability is not strictly positive, otherwise S q would diverge; such care is not necessary for 
q > 0; due to this property, the entropy is said to be expansible for q > 0). This expression 
recovers the usual BG entropy (— kY^=\P% In Pi) in the limit q — > 1. The entropic index 
q (intimately related to and determined by the microscopic dynamics, as we shall mention 
later on) characterizes the degree of nonextensivity reflected in the following pseudo-additivity 
entropy rule 

S q {A + B)/k = [S q (A)/k] + [S q {B)/k] 

+ (l-q)[S q (A)/k)[S q (B)/k], (2) 

where A and B are two independent systems in the sense that the probabilities of A + B 
factorize into those of A and of B (i.e., Pij(A + B) = pi(A)pj(B)). We immediately see 
that, since in all cases S q > [nonnegativity property), q < 1, q = 1 and q > 1 respectively 
correspond to superadditivity (superextensivity), additivity (extensivity) and subadditivity 
(subextensivity). Eq. (2) exhibits a property which has apparently never been focused 
before, and which we shall from now on refer to as the composability property. It concerns 
the nontrivial fact that the entropy S(A + B) of a system composed of two independent 
subsystems A and B can be calculated from the entropies S(A) and S(B) of the subsystems, 
without any need of microscopic knowledge about A and B, other than the knowledge of 
some generic universality class, herein the nonextensive universality class, represented by 
the entropic index q, i.e., without any knowledge about the microscopic possibilities of A 
and B nor their associated probabilities. This property is so obvious for the BG entropic 
form that the (false!) idea that all entropic forms automatically satisfy it could easily install 
itself in the mind of most physicists. To show counterexamples, it is enough to check that 
the recently introduced Anteneodo- Plastino's |1| and Curado's |17| entropic forms satisfy 



a variety of interesting properties, and nevertheless are not composable. 

The above pseudo-extensivity property can be equivalently written as follows: 

ln[l + (1 - q)S q (A + B)/k] _ ln[l + (1 - q)S q (A)/k] , ln[l + (1 - q)S q (B)/k] 



(3) 



We come back onto this form later on in connection with Renyi's entropy. 

Another important (since it eloquently exhibits the surprising effects of nonextensivity) 
property is the following. Suppose that the set of W possibilities is arbitrarily separated 
into two subsets having respectively Wl and Wm possibilities (Wl + Wm = W). We define 
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Pl = Ya=\ Pi and p M 
established that 



Y?i=w L +iPii hence pl + Pm = 1- It can then be straightforwardly 



S. 



<({Pi}) = S q (p L ,p M ) +Pl S q({Pi/PL}) + PM S q ({pi/p M }) 



(4) 



where the sets {pi/pi} and {Pi/pm} are the conditional probabilities. This would precisely 
be the famous Shannon's property were it not for the fact that, in front of the entropies 
associated with the conditional probabilities, appear p\ and p q M instead of pl and pu ■ This 
fact will play, as we shall see later on, a central role in the whole generalization of ther- 
mostatistics. Indeed, since the probabilities {p{\ are generically numbers between zero and 
unity, p\ > pi for q < 1 and p\ < pi for q > 1, hence q < 1 and q > 1 will respectively 
privilegiate the rare and the frequent events. This simple property lies at the heart of the 
whole proposal. Santos has recently shown |18[ , strictly following along the lines of Shannon 
himself, that, if we assume (i) continuity (in the {pi}) of the entropy, (ii) increasing mono- 
tonicity of the entropy as a function of W in the case of equiprobability, (iii) property (2), 
and (iv) property (4), then only one entropic form exists, namely that given in definition 
(1). Of course, the generalization of Eq. (4) to the case where, instead of two, we have 
R nonintersecting subsets {W\ + W 2 + ... + Wr = W) is straightforward fL9f . To be more 
specific, if we define 



H Pi (i = 1 > 2 , •••,#) 

Wi terms 



(5) 



(hence SjLi flj — 1)) Eq. (4) is generalized into 



R 



3=1 



(6) 



where we notice, in the last term, the emergence of what we shall soon introduce generically 
as the unnormalized g-expectation value (of the conditional entropies S q ({pi/7ij}), in the 
present case). 

Another interesting property is the following. The Boltzmann-Gibbs entropy Si satisfies 
the following relation: 



- k 



d w 

— Yp a 



-kJ^Pi ln Pi = Si 



(7) 



a=l 



i=l 



Moreover, Jackson introduced in 1909 [^(J the following generalized differential operator 
(applied to an arbitrary function f(x)): 



D q f{x) 



fjqx) - f(x) 

qx — x 



(8) 



which satisfies Di = lim q ^iD q = Abe []2T| recently remarked that 



1 _ vw rf 
, 1 2^j=i Pi 

q-1 



- k 



i=i 



a=l 



= s n 



(9) 
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This property provides some insight into the generalized entropic form S q . Indeed, the 
inspiration for its use in order to generalize the usual thermal statistics came from mul- 
tifractals, and its applications concern, in one way or another, systems which exhibit scale 
invariance. Therefore, its connection with Jackson's differential operator appears to be 
rather natural. Indeed, this operator "tests" the function f(x) under dilatation of x, in 
contrast to the usual derivative, which "tests" it under translation of x. 

Another property which no doubt must be mentioned in the present introduction is 
that S q is consistent with Laplace's maximum ignorance principle, i.e., it is extremum at 
equiprobability (p^ = 1/W Vz). This extremum is given by 

W x ~ q - 1 

S q = k (W > 1) (10) 

which, in the limit q — ► 1, reproduces Boltzmann's celebrated formula S = kin W (carved 
on his marble grave in the Central Cemetery of Vienna). In the limit W — > oo, S g diverges 
if q < 1, and saturates at k/(q — 1) if q > 1. 

Finally, let us close the present set of properties by reminding that S q has, with regard to 
{pi}, a definite concavity for all values of q (S q is always concave for q > and always convex 
for q < 0). In this sense, it contrasts with Renyi's entropy S q = (In S^ 1 pf)/(1 — q) — 
{In [1 + (1 — q)S q /k]}/(l — q), which does not have this property for all values of q . 

Before addressing other relevant quantities, let us introduce the following convenient 



functions [221: 



e^[l + (l-g)x] 1 /( 1 ^, W(x,q) (11) 

(hence, ef = e x ) with the definition supplement, for q < 1, that e x q = if 1 + (1 — q) x < 0, 
(and analogously, for q > 1, e x q diverges at x = l/(q — 1)) and 

In, x= [x^-iyil-q], V(x,g) (12) 

(hence, ln x x = In x). We can verify easily that 

e^ x =\n q e*=x, V(x,g). (13) 

For instance, Eq. (10) can be rewritten in the following Boltzmann-like form: 

S q = k \n q W (14) 

Let us also introduce the following unnormalized q-expectation value: 

w 

(A) 9 = 5>?A (15) 

i=l 

hence (A)i corresponds to the standard mean value of a physical quantity A. 

If our system is a generic quantum one, its probabilistic description is given by the density 
operator p, whose eigenvalues are the {pi}- Then, the generalized entropy is given by 

I — Tr o q 

S q = k l 1TP (Trp = l) (16) 
q — 1 
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and the unnormalized g-expectation value of an observable A which does not necessarily 
commute with p is given by 

(A) q = Tr p«A (17) 

Eq. (16) can be rewritten as follows: 

S q = -k (\n q p) q (18) 

and also as 

S q = -k (ln 2 _ 9 p) x (19) 

If our system is a generic classical one, the relevant variables are typically continuous 
variables, and its probabilistic description is given by a distribution of probabilities p(r), 
where r is a dimensionless variable in say a many-body phase space. Then, the generalized 
entropy is given by 

s , = t W*[p(,)]. { j dip{T) = 1) m 

and the unnormalized g-expectation value of an observable A(r) is given by 

(A) q = fdr [p(r)Mr) (21) 

Although we shall, in what follows, be illustrating the present formalism with the case of W 
discrete microscopic possibilities, the generic quantum and classical discussions follow along 
the same lines, mutatis mutandis. 



B. CANONICAL ENSEMBLE 



Once we have a generalized entropic form, say that given in Eq. (1) (or an even more 
general one, or a different one), we can use it in a variety of ways. For instance, if we are inter- 
ested in information theory, some optimization algorithms, image processing, among others, 
we can take advantage of a particular form in different ways. See, for instance, |TT|fiTS| , |2"B| , |2~3 



and references therein, where it can be verified that not less than 25 (!) different entropic 
forms have received, along the years, a great variety of technological and mathematical ap- 
plications. For instance, the Renyi entropy mentioned above has been quite useful in the 
geometrical characterization of strange attractors and similar multifractal structures (see 



25|| and references therein). 

However, if our primary interest is Physics, this is to say the (qualitative and quanti- 
tative) description and possible understanding of phenomena occurring in Nature, then we 
are naturally led to use the available generalized entropy in order to generalize statistical 
mechanics itself and, if unavoidable, even thermodynamics. It is along this line that we 
shall proceed from now on (see also p6| ). To do so, the first nontrivial (and quite ubiqui- 
tous) physical situation is that in which a given system is in contact with a thermostat at 
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temperature T. To study this, we shall follow along Gibbs' path and focus the so called 
canonical ensemble. More precisely, to obtain the thermal equilibrium distribution associ- 
ated with a conservative (Hamiltonian) physical system in contact with the thermostat we 
shall extremize S q under appropriate constraints. These constraints are []i~5 



and 



i=l 



(norm constraint) 



yw n i f . 



T w n q 
2^i=i Pi 



U q (energy constraint) 



(22) 



(23) 



where {e^} are the eigenvalues of the Hamiltonian of the system. We shall refer to ((...)) q 
as the normalized q- expectation value and to U q as the generalized internal energy (assumed 
finite and fixed). It is clear that, in the g — ^ 1 limit, these quantities recover the standard 
mean value and internal energy respectively. We immediately verify that, for any observable, 



«■■•»< 



(1), 



(24) 



The outcome of this optimization procedure is given by 



Pi 



l-(l- q )f3(e i -U q )lT,U (P. 



1-9 



(25) 



with 



w 



zM = E 



w 



1 - (1 - qtfte - U q )f£ ( Pj 



3=1 



1-9 



(26) 



It can be shown that, for the case q < 1, the expression of the equilibrium distribution is 
complemented by the auxiliary condition that Pi = whenever the argument of the function 



becomes negative (cut-off condition). Also, it can be shown |L5| that 



l/T = dS q /dU q , Wq (T = l/(kP)). 



(27) 



Furthermore, it is important to notice that, if we add a constant eo to all {e^}, we have (as it 
can be self-consistently proved) that U q becomes U q +e , which leaves invariant the differences 
{e.j — U q }, which, in turn, (self-consistently) leaves invariant the set of probabilities {pi}, 
hence all the thermostatistical quantities. It is also trivial to show that, for the independent 
systems A and B mentioned previously, U q (A + B) = U q (A) + U q (B), thus recovering the 
same form of the standard (q = 1) thermodynamics. 
It can be shown that the following relations hold: 



Efe) 9 = (z q 



(28) 



i=l 



7 



(29) 



and 



a (z Y~i - i 



where 

(z q y-"-i _ (z q y-"-i 

1 — q 1 — q 



PU q . (31) 



Let us now make an important remark. If we take out as factors, in both numerator and 
denominator of Eq. (25), the quantity 1 + (1 — q)PU q /Y^fLi (Pj) 9 , and then cancel them, 
we obtain 



[l-a-g)/^] — 



ir 



with 



pAP) = 1 - — v z r " K^EIi- (i-^l 1 ^ I 02) 



^ EW + V«)/*/. <r ^ 1/(W) <33) 



where /?' is an increasing function of (3 [271 . 

Let us now address the all important question of the connection between experimental 
numbers (those provided by measurements), and the quantities that appear in the theory. 
The definition of U q suggests the following normalized q- expectation values 



O q = m)) q = ^ Pl 7 (34) 
2^i=i Pi 



where O is any observable which commutes with the Hamiltonian, hence with p. If it does 
not commute, Eq. (34) is generalized into 

Consistently, O q is the mathematical object to be identified with the numerical value pro- 
vided by the experimental measure. Later on, we come back onto this crucial point. 



At this point let us make some observations about the set of escort probabilities |28 
{Pi^} defined through 

q W 

(£^ = 1) (36) 
Z^=i Pj i=i 

from which follows the dual relation 
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The W = 2 illustration of is shown in Fig. 1. As anticipated, q < 1 (q > 1) privileges 
the rare (frequent) events. 

Let us first comment that Eqs. (36) and (37) have, within the present formalism, a role 
somehow analogous to the direct and inverse Lorentz transformations in Special Relativity 
(see p9[ and references therein). Second, we notice that O q becomes an usual mean value 
when expressed in terms of the probabilities {P^}, i. e., 



0^%^ = E^- (38) 



and 



w 

^P? ) e i = U q . (39) 



i=l 

The final equilibrium distribution reads 

F(« = P -(!-«>/*■]*, . (40 ) 

Er=i[i-(i-9).3'«]~ 

If the energy spectrum {e^} is associated with the set of degeneracies {gi}, then the above 
probability leads to the following one (associated with the level and not the state i) 

p(«) = g .n-(i-3)/^ ^ (41) 

T,all levels 9k[l ~ (1 ~ g)0' e fc] ^ 

If the energy spectrum {e^} is so dense that can practically be considered as a continuum, 
then the discrete degeneracies yield the function density of states g(e), hence 

P M=gW II- d-^l^ , (42 ) 

The density of states is of course to be calculated for every specific Hamiltonian (given the 
boundary conditions). For instance, for a c/-dimensional ideal gas of particles or quasipar- 
ticles, it is given [[3(J by g(e) oc e* _1 , where 5 is the exponent characterizing the energy 
spectrum e oc K s where K is the wavevector (e.g., 5=1 corresponds to the harmonic oscil- 
lator, 5 = 2 corresponds to a nonrelativistic particle in an infinitely high square well, etc). 
In Figs. 2 and 3 we see typical energy distributions for the particular case of a constant 
density of states. Of course, the q = 1 case reproduces the celebrated Boltzmann factor. 
Notice the cut-off for q < 1 and the long algebraic tail for q > 1. 

All the above considerations refer, strictly speaking, to thermodynamic equilibrium. The 
word thermodynamic makes allusion to "very large" (N — > oo, where N is the number of 
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microscopic particles of the physical system). The word equilibrium makes allusion to asymp- 
totically large times (t — * oo limit) (assuming a stationary state is eventually achieved). The 
question arises: which of them first? Indeed, although both possibilities clearly deserve the 
denomination "thermodynamic equiulibrium" , nonuniform convergences might be involved 
in such a way that lim^^co limt^oo could differ from limt^oo UmN^oo- To illustrate this 
situation, let us imagine a classical Hamiltonian system including two-body interactions 
decaying at long distances as l/r a in a <i-dimensional space, with a > 0. If at > d the 
interactions are essentially short- ranged, the two limits just mentioned are basically inter- 
changeable, and the prescriptions of standard statistical mechanics and thermodynamics are 
valid, thus yielding finite values for all the physically relevant quantities. In particular, the 
Boltzmann factor certainly describes reality, as very well known. But, if < a < d, nonex- 
tensivity is expected to emerge, the order of the above limits becomes important because of 
nonuniform convergences, and the situation is certainly expected to be more subtle. More 
precisely, a crossover (between q ^ 1 and q = 1 behaviors) is expected to occur at t = t(N). 
If lim^^ 00 T{N) = oo, then we would indeed have two (or even more) different and equally 
legitimate states of thermodynamic equilibrium, instead of the familiar unique state. The 
conjecture is illustrated in Fig. 4. 

A wealth of works has shown that the above described nonextensive statistical mechan- 
ics retains much of the formal structure of the standard theory. Indeed, many important 
properties have been shown to be q-invariant. Among them, it is mandatory to mention 

(i) the Legendre transformations structure of thermodynamics |14] , [15| ; 

(ii) the H-theorem (macroscopic time irreversibility), more precisely, that, in the presence 
of some irreversible physical evolution, dS q /dt > 0, = and < if q > 0, =0 and < 0, 
respectively, the equalities holding for equilibrium [|3Ty32f ; 

(hi) the Ehrenfest theorem (correspondence principle between classical and quantum me- 
chanics) f33fl ; 

(iv) the Onsager reciprocity theorem (microscopic time reversibility) pi,|55|; 

(v) the Kramers and Wannier relations (causality) ||35|| ; 

(vi) the factorization of the likelihood function (Einstein' 1910 reversal of Boltzmann's for- 
mula) 0; 

(vii) the Bogolyubov inequality [j36fl ; 

(viii) thermodynamic stability (i.e., a definite sign for the specific heat) ||37|| ; 

(ix) the Pesin equality . 

In contrast with the above quantities and properties, which are q- invariant, some others 
do depend on q, such as 

E3; 



(i) the specific heat 

(ii) the magnetic susceptibility (40 



(iii) the fluctuation-dissipation theorem (of which the two previous properties can be con- 
sidered as particular cases) f40|j ; 

(iv) the Chapman- Enskog expansion, the Navier-Stokes equations and related transport co- 
efficients PT|j ; 

(v) the Vlasov equation [|I^ , [l3|| ; 



vi) the Langevin, Fokker-Planck and Lindblad equations p4}44^]; 
stochastic resonance 



vn 



49 



(viii) the mutual information or Kullback-Leibler entropy P2|,|5U 
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A remark is necessary with regard to both sets just mentioned. Indeed, these properties 
have in fact been studied, whenever applicable, within unnormalized g-expectation values 
for the constraints, rather than within the normalized ones that we are using herein. Nev- 
ertheless, they still hold because they have been established for fixed /3, which, through Eq. 
(33), implies fixed (3' . 

Finally, let us mention various important theoretical tools which enable the thermosta- 
tistical discussion of complex nonextensive systems, and which are now available (within the 
unnormalized and/or normalized versions for the g-expectation values) for arbitrary q. We 
refer to 

i) Linear response theory ||35|| ; 

ii) Perturbation expansion |TjJ ; 

iii) Variational method (based on the Bogoliubov inequality) 



iv) Many-body Green functions |52 



v) Path integral and Bloch equation |[53|| , as well as related properties 

vi) Quantum statistics and those associated with the Gentile and the Haldane exclusion 
statistics |55|-p7|; 

(vii) Simulated annealing and related optimization, Monte Carlo and Molecular dynamics 
techniques pSHSB . 



III. THEORETICAL EVIDENCES AND CONNECTIONS 

A. Levy-type anomalous diffusion 

An enormous amount of phenomena in Nature follow the Gaussian distribution: measure- 
ment error distributions, height and weight distributions in biological individuals of given 
species, Brownian motion of particles in fluids, Maxwell-Boltzmann distribution of particle 
velocities in a variety of systems, noise distribution in uncountable electronic devices, energy 
fluctuations at thermal equilibrium of many systems, to only mention a few. Why is it so? 
Or, equivalently, what is their (thermo) statistical foundation? This fundamental problem 
has already been addressed, particularly by Montroll, and satisfactorily answered (see |7j 
and references therein) . The answer basically relies onto two pillars, namely the BG entropy 
and the standard central limit theorem. However, the Gaussian is not the only ubiquitous 
distribution: we also similarly observe Levy distributions (in micelles [^J, supercooled laser 



|70| , fluid motion |71[, wandering albatrosses [f72[ , heart beating |73[, financial data |74|-|76|, 
among many others). So, once again, what is the (thermo) statistical foundation of their 
ubiquity? This relevant question has also been addressed, once again by Montroll and col- 
laborators [0] among others. In this case however, a satisfactory answer has been missing 
for a long time. The first successful step toward (what we believe to be) the solution was 
performed in 1994 by Alemany and Zanette |77j], who showed that the generalized entropic 
form S q was able to provide a power-law (instead of the exponential-lew associated with 
Gaussians) decrease at long distances. Many other works followed along the same lines 
|78| , [79|| . In |79| it was exhibited how the Levy-Gnedenko central limit theorem also plays 
a crucial role by transforming, through successive iterations of the jumps, the power- law 
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obtained from optimization of S q into the specific power-law appearing in Levy distribu- 
tions. Summarizing, in complete analogy with the above mentioned Gaussian case (and 
which is recovered in the more powerful present formalism as the q = 1 particular case), the 
answer once again relies onto two pillars, which now are the generalized entropy S q and the 
Levy-Gnedenko central limit theorem. 



The arguments have been very recently re-worked |80| out on the basis of the normalized 
q- expectation values introduced in []I5| |. These are the results that we briefly recall here. 



Let us write S q as follows: 



1 " I-oo f W Pi*)}" 



S q \p(x)] = k 



(43) 



where x is the distance of one jump, and a > is the characteristic length of the problem. 
We optimize (maximize if q > 0, and minimize if q < 0) S q with the norm constraint 
J^oo dx p(x) = 1, as well as with the constraint 



fZ dx\p(x)]i 



a' 



(44) 



We straightforwardly obtain the following one-jump distribution. 
If q > 1: 



If q = 1: 



If q < 1: 



Pg(X) 



a 



q-1 



7T 



1/2 



i 



• 2 1 V(9-l) 



(45) 



p q (x) 



1 r 1 
a 



2tt 



1/2 



-{x/of/2 



(46) 



Pg(X) 



a 



1-q 



7T 



(3-g) 



V2 r (ji^)) 

r(fEf) 



1-? r 
3 — g a 2 



V(l-9) 



(47) 



if |x| <a[(3-g)/(l-g)] 1 / 2 and zero otherwise. 

We see that the support of p q (x) is compact if q G (— oo,l), an exponential behavior 
is obtained if q = 1, and a power-law tail is obtained if q > 1 (with p g (x) oc (cr/x) 2 /^ 1 ) 
in the limit |x|/cr — > oo). Also, we can check that ((x 2 ))i = {x 2 )i = J^dx x 2 p q (x) is 
finite if q < 5/3 and diverges if 5/3 < q < 3 (the norm constraint cannot be satisfied if 
q > 3). Finally, let us mention that the Gaussian (q = 1) solution is recovered in both limits 
q — >• 1 + and q — > 1 — by using the q > 1 and the q < 1 solutions respectively. This 
family of solutions is illustrated in Fig. 5. 
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We focus now the iV-jump distribution p q (x,N) = p q (x) * p q (x) * ... * p g (x) (iV-folded 
convolution product). If q < 5/3, the standard central limit theorem applies, hence, in the 
limit N — > oo, we have 



Pq(x, N) 



1 



a 



5-3q 



2tt(3 - q)N 



1/2 



exp 



5 — 3q x 1 
' 2(3 -q)No* 



(48) 



i.e., the attractor in the distribution space is a Gaussian, consequently we have normal 
diffusion. If, however, q > 5/3, then what applies is the Levy-Gnedenko central limit 
theorem, hence, in the limit N — > oo, we have 



Pq{N,x) ~ L 7 (x/N 



where L 7 is the Levy distribution with index 7 < 2 given by 



7 



3-q 
q-1 



(5/3 < q < 3) 



(49) 



(50) 



Through the Fourier transforms of both Eq. (48) and (49), we can characterize the width 
A q (dimensionless diffusion coefficient) of p g (x, N). We obtain 



5-3g 



(q < 5/3) 



and 



2 


r«- 1- 


3-q 

2(5-1) y 


r 3q - 5 -I 


Tpl/2 


-3-q- 




-2(q-l)\ 



(5/3 < q < 3) 



(51) 



(52) 



These results are depicted in Fig. 6. This result should be measurable in specifically devised 
experiments. More details can be found in [ 30 1 and references therein. What we wish to 
retain in this short review is that the present formalism is capable of (thermo) statistically 
founding, in an unified and simple manner, both Gaussian and Levy behaviors, very ubiqui- 
tous in Nature (respectively associated with normal diffusion and a certain type of anomalous 
superdiffusion). 



B. Correlated-type anomalous diffusion 

There are some phenomena exhibiting anomalous (super and sub) diffusion of a type 
which differs from the one discussed in the previous subsection. We refer to the so called 
correlated-type of diffusion. We consider here a quite large class of them, namely those 
associated with the following generalized, Fokker-Planck-like equation: 

I cp(*, t)r = - JVw frfo + d^hx, or (53) 

where (/i, v) G TZ 2 , D is a dimensionless diffusion-like constant, F(x) = —dV/dx is a di- 
mensionless external force (drift) associated with a potential V(x), and (x,t) is a dimen- 
sionless 1 + 1 space-time. If p, — 1, we can interpret p(x,t) as a probability distribution 
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since / dx p(x,t) = 1, Vt can be satisfied. If 7^ 1, then p(x,t) must be seen as a 
density function. The word "correlated" is frequently used in this context due to the fact 
that D(d 2 /dx 2 )\p(x, t)] u = (d/dx){Du[p(x,t)] u ~ 1 (d/dx) p(x,t)}, i.e., an effective diffusion 
emerges, for v ^ 1, which depends on p(x, t) itself, a feature which is natural in the presence 
of correlations. The ji — 1 particular case of this nonlinear equation is commonly denomi- 
nated "Porous medium equation", and corresponds to a variety of physical situations (see 
4"5|] and references therein for several examples). 

The first connection of Eq. (53) with the present nonextensive statistical mechanics was 



established in 1995 by Plastino and Plastino |45|| . They considered a particular case, namely 
fM — 1 and F(x) = —k 2 x with k 2 > (so called Uhlenbeck-Ornstein processes), and found 



an exact solution which has the form of Eq. (43-45). Their work was generalized in f46 
where arbitrary /i and F(x) = ki — k 2 x were considered. The explicit exact solution of 
Eq. (53), for all values of (x,t), was once again found by proposing an Ansatz of the form 
of Eqs. (45-47), i.e., the form which optimizes S q with the associated simple constraints. 
This form eventually turns out to be the Barenblatt one, useful in related problems. Here, 
let us restrict ourselves to just reproduce the exact solution of Eq. (53) assuming that 
p(x, 0) = S(x), this is to say, a Dirac delta distribution. We obtain |46| 



^^ilzMyig (54) 

where 

q — 1 + [A — v (55) 



and 



m l z q (t) 

with 



(56) 



Z q {t)=ZM[{l-±y' r + ± 1 '™ , (57) 



K2 ~ 2i/D/3(0)[Z,(0)]"-" (58) 

and 

T = 7^-^ (59) 

Summarizing, by using the form which optimizes S g , it has been possible to find the phys- 
ically relevant solution of a nonlinear equation in partial derivatives with integer derivatives. 



It can be shown |8T[ that the problem that was solved in the previous subsection corresponds 
to a linear equation in partial derivatives but with fractional derivatives. We believe that we 
are allowed to say that an unusual mathematical versatility has been observed, within the 
present nonextensive formalism, in this couple of nontrivial examples of anomalous diffusion. 
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C. Stellar polytropes and other self-gravitating systems 



The present formalism has been applied to a variety of astrophysical [^,0 and cosmo- 



logical |82[ self-gravitating systems. In some sense, this is something natural to do given the 



long range of the gravitational interaction. This was, in fact, the first physical application 
of nonextensive statistics. We do not intend here to reproduce details. Our present aim is to 
remind that it is well known in astrophysics that, within the standard thermodynamical ap- 
proaches, it is not possible to simultaneously have finite values for the total energy, entropy 
and mass of a self-gravitating system. Plastino and Plastino were the first to show, in 1993, 
that this physically desirable situation can be achieved if we allow q to sufficiently differ 
from unity ! In fact, it can be shown (by considering the Vlasov equation in .D-dimensional 
Schuster spheres) that the problem becomes a mathematically well posed one if q < q*, 
where the critical value q* is given by 

8-(D-2) 2 

Q ~ 8- (D- 2)2 + 2(^-2) ' (60) 

For D = 3 we recover the 7/9 relatively known value. Also, we notice that D = 2 implies 
q* = 1, which is very satisfactory since it is known that D < 2 gravitation is tractable within 
standard thermodynamics. 



D. Zipf-Mandelbrot law 

The problem we focus here first appeared in Linguistics. However, its relevance is quite 
broad, as it will soon become clear. Suppose we take a given text, say Cervante's Don 
Quijote, and order all of its words from the most to the less frequent; we refer to the ordered 
position of a given word as its rank R (low rank means high frequency uj of appearance in 
the text, and high rank means low frequency). Zipf |||] discovered that, in this as well as 
in a variety of similar problems, the following law is satisfied: 



A Br* (Zipf law) (61) 



where A > and £ > are constants. Later on, Mandelbrot |73| suggested that such 
behavior was reflecting a kind of fractality hidden in the problem; moreover, he suggested 
how the Zipf law could be numerically improved: 

u = — ^ {Zipf — Mandelbrot law) (62) 

This expression has been useful in a variety of analysis. The connection we wish to mention 



here is that in 1997 Denisov [84] showed that, by extending (to arbitrary q) the well known 
Sinai-Bowen- Ruelle thermodynamical formalism of symbolic dynamics (i.e., by considering 
S q instead of Si), the Zipf-Mandelbrot law can be deduced. He obtained 
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and D = d/(q — 1) (d being a positive constant), i.e., 



u a [i + 5 - 1) (g>1) (64) 

Clearly, to make the discussion complete, a model would be welcome, which would provide 
quantities such as q and d. Nevertheless, Denisov's arguments have the deep interest of 
explicitly exhibiting that the Zipf-Mandelbrot law can be seen as having a nonextensive 
foundation. Fittings with experimental data will be shown later on in connection with the 
citations of scientific papers. 



E. Theory of financial decisions; Risk aversion 

An important problem in the theory of financial decisions is how to take into account 
extremely relevant phenomena such as the risk aversion human beings (hence financial op- 
erators) quite frequently feel. This kind of problem has, since long, been extensively studied 
by Tversky |85j and co-workers. The situation can be illustrated as follows. What do you 



prefer, to earn 85,000 dollars or to play a game in which you have 0.15 probability of earning 
nothing and 0.85 probability of earning 100,000 dollars ? In fact, most people prefer take the 
money. The problem of course is the fact that the expectation value for the gain is one and 
the same (more precisely 85,000 dollars) for both choices, and therefore this mathematical 
tool does not reflect reality ! The same problem appears if one expects to loose 85,000 and 
the chance is given for playing a game in which, if you win, you pay nothing, but, if you 
loose, you pay 100,000 dollars. In this case, most people choose to play. So, the experimental 
facts are that most human beings are risk-averse when they expect to gain, and risk-seeking 
when they expect to loose ! The problem is how to put this into mathematics. 



values as follows: 



Following |86||, let us introduce, for the above gain problem, normalized g-expectation 

({gain))f kethemoney = 85,000 (65) 



and 

(( . ))P lay the game = 100, 000 X 0.85" + X 0. 15" = 100,000 x 0.85" 

Uy ' 1q 0.859 + 0.159 0.859 + 0.159 1 ' 

Since most people would prefer the money, this means that most people have q < 1 for this 
particular decision problem. 
For the loss probem we have: 

((gain))f kethemoney = -85,000 (67) 

and 

(( . ))P lay the game = -100,000 x 0.859 + x 0.159 = -100,000 x 0.859 

[[g ))q 0.859 + 0.159 0.859 + 0.159 1 ' 

Since in this case most people would prefer to play, this means that, consistently with the 
previous result, most people have q < 1 for the particular decision problem we are considering 
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now. In some sense, we have some epistemological progress ! Indeed, the statement "most 
people have (for this type of amount of money) q < 1", unifies the previous two separate 
statements concerning expectation to gain and expectation to loose. 

Let us address now the following question: how can we measure the value of q associated 
with a particular individual ? We illustrate this interesting point with the example of the 
gain. The person is asked to choose between having V dollars or playing a game in which, 
if the person wins, the prize will be 100, 000 dollars and, if the person looses, he (she) will 
receive nothing. As before, the person is informed that his (her) probability of winning is 
0.85 (hence, the probability of loosing is 0.15). Then we keep gradually changing the value 
V and asking what is the preference. At a certain critical value, noted V c , the person will 
change his (her) mind. Then, the value of q to be associated with that person, for that 
problem, is given by the following equality 

100,000 x 0.85 9 

' = Vc (69 

0.85^ + 0.159 v ; 

(See Fig. 7). The ideally rational operator corresponds to q — 1. For this gain problem, 
the risk-averse operators correspond to q < 1, and the risk-seeking ones to q > 1. It is clear 
that models for stock exchange can be formulated by using these remarks. Such an effort is 



presently in progress [37 . 



F. Physiology of vision 

Physiological perceptions such as the visual perception are since long known to focus upon 



rare events (e.g., a red spot on a white wall). Barlow |p8 |, among others, has recurrently 
stressed our attention on the fact that, at the action decision level, the various possibilities 
should enter with a weight proportional to — lnpj, and not proportional to pi, pi being the 
a priori probability of occurrence of that particular event; indeed, — lnp« diverges when 
Pi — > 0. He even argues that evolutionary arguments hold very well together with such 
hypothesis. To privilege rare events is precisely what happens, in the present formalism, 



whenever q < 1. Let us be more specific: if we consider the < q « 1 limit, we obtain [59 



<1 W 



1 VV H 

S q /k = ~^ Pi ~ W - 1 + q[W - 1 - £(- Inft)], (70) 
1 1 ?;=i 



W WW 



(0) q = J2p!O^J20i- q E(- ^Pi) Oi (71) 



i=l i=l i=l 



and 



«°»« = ^l 1 + 4 w yZ.o, J / (72) 

where O is an arbitrary observable. Leaving aside several constant quantities that appear 
above, we immediately observe the prominent role which — In pi plays in these expressions. 
Consistently, the g — ^ limit of the present formalism could well be of some utility in the 
theoretical analysis of the physiological phenomena focused here. 
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IV. EXPERIMENTAL EVIDENCES AND CONNECTIONS 



A. D = 2 turbulence in pure-electron plasma 

A few years ago, in 1994, Huang and Driscoll |j) exhibited some quite interesting non- 
neutral plasma experimental results obtained in pure-electron plasma (confined in a 20 cm 
long and 6 cm wide metallic cylindric Penning trap with a 10~ 10 torr vacuum in its interior) 
in the presence of an external axial magnetic field (507 Gauss). In the interval 2-100 ms 
after every single electric shot (generating the electron plasma), it was observed a turbulent 
axisymmetric metaequilibrium state, the electronic density radial distribution of which was 
measured. Its average (over typically 100 shots) monotonically decreased with the radial 
distance, disappearing at some radius sensibly smaller than the radius of the container (i.e., 
a cut-off was observed). The experiment was recently redone |90| under slightly modified ex- 



perimental conditions (a slow external rotation was imposed in such a way as to compensate 
the small energy dissipation existing in the plasma), and essentially the same metaequi- 
librium state was observed during lapses of time as long as 27 hours, or even longer ! In 
addition to the 1994 experiment, the authors also proposed M a phenomenological theory 
trying to reproduce the experimentally observed profile. Their proposal consisted on the op- 
timization, for a given model, of a functional of the electron density p(r) under constraints, 
namely conservation of total mass, angular momentum and energy. They presented four 
different attempts. The first one (Point Vortex Maximum Entropy) consisted in optimizing, 
for a point vortex representation of the plasma, the BG entropy: it failed in reproducing the 
experimental data. The second attempt (Fluid Maximum Entropy) was essentially the same 
as the previous one, but using a fluid model for the plasma: the failure was even bigger. 
They assumed next that the problem possibly relied, not so much in the particular plasma 
model, but rather in the chosen functional to be optimized. In their third attempt (Global 
Minimum Enstrophy) , they turned back to the point vortex model, but optimized the enstro- 
phy instead of the BG entropy. The result was better than the two first attempts, but had 
the physically unacceptable feature of producing a negative electron density at sufficiently 
high radius. They then addressed their fourth attempt (Restricted Minimum Enstrophy), 
whose only difference with the third one was the fact of introducing an out-of-the pocket 
cut-off of the electron density at the proper value of the radius. This procedure was, finally, 
successful, and a very good first-approximation fitting was obtained ! The effort done by 
Huang and Driscoll was, on top of the high merit of a remarkable experiment, extremely 
pedagogical and elucidating: the main theoretical problem was not the model, but rather 
the choice of the functional to be optimized, i.e., the statistics. 

The next important step in this story was done by Boghosian. He realized in 1995 
and published JEJ in 1996 that the Huang and Driscoll fourth, successful attempt precisely 
corresponds to the optimization of S q with q = 1/2 ! Indeed, by following the recipes of the 
present generalized thermostatistics, he re-obtained, for the electron density profile, the same 
differential equation produced within the Restricted Minimum Enstrophy phenomenological 
theory, with the supplementary bonus of not having to introduce in an ad hoc manner the 
necessary cut-off. Indeed, as already argued, all q < 1 cases exhibit a cut-off intrinsic to the 
formalism, and the radial position of that cut-off nicely fits the experimental value. 
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The next step was performed in 1997 by Anteneodo and myself |9T[| (in fact, after related 
remarks by Boghosian himself). The Restricted Minimum Enstrophy theory is based on 
the enstrophy functional, which belongs to the general discussion of Casimir invariants; its 
form is in fact that of the order 2 Casimir invariant. Consequently, an epistemo logically 
conservative theoretical viewpoint is to appreciate Boghosian's effort as just a formal inter- 
esting remark, with no real physical necessity. It happens, however, that, for r — > r c — 0, 
(r c = cut-off radius) the enstrophy theory yields p(r) oc (r c — r) whereas the experimental 
data fit much better a vanishing derivative at r c ! We followed |)1|] along Boghosian's lines 
and generalized his theory for arbitrary q. We obtained the generalized differential equation 
for p(r) and showed that p(r) cx (r c — r) 9 ^ 1 " 9 ). Consequently, the experimental data fit 
better for q slightly above 1/2. This, together with the numerical solution of the differential 
equation, advanced q ~ 0.55 as a better value for satisfactory overall fitting. (Better fittings 
would probably demand for a model more sophisticated than the point vortex one used 
here). The conceptually important point of this discussion is that Casimir invariants are 
characterized by integer exponents (in p(r)), hence none of them can be related to a value of 
q close to 0.55. From this standpoint, the present formalism appears as the only satisfactory 
phenomenological theoretical approach available in the literature at the present time. 

The last step of this analysis was performed very recently by Anteneodo [ 9~2"] . Indeed, the 
calculations above recalled [43|,91] were done by using unnormalized g-expectation values. 
However, as already mentioned and used in the present review, it has been recently argued 



[ 15] that normalized g-expectation values should be used instead. It is therefore important 
to check that the present discussion and results for turbulence remain essentially invariant. 
This is now done []92 



and it is this theory that we present in what follows. 



The generalized entropy and associated constraints are given by 



s q [g] 



q-l 



J(g- g q )d\ 



(73) 



/ gd 2 r = 1 (mass conservation) 

L q = L (angular momentum conservation) 

U q = U (energy conservation) , 



Jg q d 2 r 



jgid 2 r 

where g(r) is the probability distribution. Moreover, the scaled electrostatic potential 
<P(r) _ f g q (r')G(r,r')d 2 r' 



(74) 



Jg g d 2 r 



with V 2 G(r,r') =4n5(r-r'), 



satisfies 



V 2 — = 4tt 



g q 

I g q d 2 r' 



(75) 



(76) 
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The constrained optimization of S q [g] (S(S g — a J gd 2 r — XL q — (3U q )) now yields 

1 — Q9 q 1 A 2 x /3 g x A ,_i /_ x 

g -l " a " iV gr ^ + iV V 9 ^ iV - 2U N )9 « = ° (77) 

(where / g q d 2 r = N) or 

9 l q ~ q - g i_, A 2 3 <f) q A /3 

T^r " ag iv gr + ]v V + q{ n - 2 V = (78) 

or, taking the Laplacian of both sides, 

[1 + «(1 " " 4 ^ + 47r ^^ = ( 79 ) 



which can be rewritten as 



9 q -^ + 9 j = 9l(B^ q - ^) (80) 

9q ' 



where A f = 4g^/[l + a(l - g)] and 5 1 " = 4vrg^/[l + a(l - g)]. 
Alternatively, identifying p q = g q /N, we have 



P q ~ 2 -^& + i = q p^{Bp q - A), (81) 
q Pq r 



with A = A*N i and S = B'N * . This equation precisely is the one appearing in [91 



which, for q = 1/2, recovers that of f4"3|| . For any chosen q, the values of the parameters 
(A, B) are obtained by imposing the experimental values of total angular momentum and 
energy This phenomenological theory has, therefore, only one fitting parameter (q). As said 
before, q = 1/2 exactly reproduces the Huang and Driscoll's Restricted Minimum Enstrophy 
profile. The best overall fitting is, however, obtained for a value of q slightly above 1/2. 



B. Solar neutrino problem 

As easily conceivable, the core of the Sun is a very complex and turbulent plasma, 
within which an enormous amount of nuclear reactions take place. Many of them constitute 
chains of nuclear reactions in which neutrinos are produced. For instance, the p-p chain is 
described in |93| . Through a quite complete analysis of the production of neutrinos within 



the so called Solar Standard Model (SSM), it is possible to predict the neutrino flux onto the 
Earth. However, the actual flux measured in a variety of underground laboratories (Gallex, 
Sage, Kamiokande, Super-Kamiokande, Homestake) roughly amounts to only half of the 
predicted value ! This problem is currently referred to as the "solar neutrino problem". 
Two nonexclusive sources of explanation of this enigmatic discrepancy are: (i) the possible 
neutrino oscillations, which would make that only part of the predicted value would be 
detectable on the Earth; (ii) the current use of the SSM might be incorrect because it uses 
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BG thermal statistics, which could be inappropriate for the solar plasma. Clayton |L0| was 
the first to address the second possibility, as far as 25 years ago ! Indeed, he assumed an 
hypothetic distribution of energies essentially given by 

p{E) oc e-W e~W E)2 (82) 

The particular value 5 = obviously recovers BG statistics. Clayton showed that a small 
value of 8 (5 ~ 0.01) was enough to make the theory consistent with the experimental data 
that were available at that time. Quarati and co-workers remarked (preliminarily in 1996 
and in more refined calculations since then [^]) that, since the needed 5 is very small, 



the ansatz distribution could as well be the power-law one which appears in the present 
formalism. By identifying the first corrections (to BG) of both distributions, they obtained 

t = ^ (83) 

Consequently, values of q quite close to unity are enough to fit the solar neutrino discrepancy. 
Once again, we verify the extreme efficiency that modifications of the statistics can have. 



C. Peculiar velocities in Sc galaxies 



From the data obtained by the Cosmic Background Explorer (COBE), it has been pos- 
sible to infer the distribution of peculiar velocities of certain groups of spiral (Sc) galaxies 
(we recall that by peculiar velocity we mean the residual velocity after the global universe 
expansion velocity has been substracted). Bahcall and Oh O] developed four theoretical 
attempts (namely Cold Dark Matter with Q = 0.3 and with Q = 1.0, Hot Dark Matter with 
Q = 1.0 and Primeval Barionic Isotropic with Q = 0.3). All the attempts were done within 
BG statistics. The less unsatisfactory fitting was obtained for the CDM model with Q = 0.3. 
In fact, all the attempts exhibit a long tail towards high velocities, whereas the experimental 
data show a pronounced cut-off at about 500 Km s -1 . It is relevant to mention that all the 
models that were used had several fitting parameters, and nevertheless could not get rid of 
the tail. A fitting was then advanced using the present formalism with only two free 
parameters, one of them being q and the other one a characteristic velocity. The function 
that was used was the g-generalized Maxwell distribution, essentially corresponding to an 
ideal classical gas. The quality of the fitting is quite remarkable, far better than those corre- 
sponding to the already mentioned four attempts. Once again, one sees that modifications 
of the statistics can be sensibly more efficient than modifications of the model. A famous 
example along this line is provided by the completely different physics associated with a gas 
of free fermions or of free bosons, i.e., a Fermi-Dirac ideal gas or a Bose-Einstein ideal gas 
(same model but different statistics). 



D. Nonlinear inverse bremsstrahlung absorption in low pressure argon plasma 



Liu et al |12| provided in 1994 strong evidence of the existence of non-Maxwellian velocity 
distributions in a specific plasma experiment, where low pressure argon is exposed to pulsed 
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discharges. During the afterglow, measurements of the inverse bremsstrahlung of intense 
microwaves is performed. The experimental setting is such that Coulombian collisions are 
dominant. The experimental data were fitted with the following flat-topped distribution: 

f(v) ex exp[-(v/v m ) m ] (84) 

with m > 2. Souza and myself |J7J showed in 1997 that the same data can equally well be 
fitted with 

f(v)cx[l-(l- q )(v/v g )Y {1 - q) (85) 

with q > 1. Furthermore, if we expand both fitting functions in the neighborhood of the 
Gaussian case, we obtain that 

Q = -J ( 86 ) 

In both fittings, the exponents m and q depend on the microwave power. In order to 
discriminate between the two fitting functions, quite precise and systematic experiments 
would be needed, in particular exploring the actual dependence of the results on the power. 



E. Cosmic background microwave radiation 

The most accurate data concerning the cosmic microwave background radiation have 
been obtained with the FIRAS (Far-infrared absolute spectrophotometer) instrument in the 
COBE (Cosmic background explorer) satellite |9l| . These data are known to follow, in the 
2 — 20 cm" 1 region, Planck's black-body law. In 1995, Sa Barreto, Loh and myself ||99|| , 



as well as Plastino, Plastino and Vucetich ||100|| (and several others since then), analyzed 
within what precision one is allowed to assume q = 1. The result that has systematically 
come out from these analyses is \q — 1| < 10" 4 . If new observations were performed in the 
future which would be say 10 times more precise than the available ones [ [9^1 , this bound 



would be attained. Consequently, we would know better within what degree of confidence 
extensive thermostatistics can be used for this cosmological problem. If q ^ 1 turns out to 
be clearly confirmed, it is not excluded that we would have to revise our notions about the 
structure of space-time at the appropriate scales (possibly, Planck's length). It might come 
out that the physics at that level are better described by /imte-difference equations than by 
differential equations ! 



F. Electron-positron collisions 

The electron-positron annihilation into a virtual photon and the subsequent creation of 
a quark-antiquark pair provides the cleanest environment for the hadroproduction. Each 
of the two initial partons begins a complex cascade related to the strong-coupling long- 
distance regime of Quantum Chromodynamics. A partially successful global description of 
the hadroproduction has been provided through a thermodynamical equilibrium approach, 
mainly that of Hagedorn in 1965 [ |101|| . This theory provides the following prediction: 
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1 da 3/2 



cpt' exp -p T /T {p T > T ) (87) 



where a is the distribution of the transverse momenta px, To is a characteristic temperature 
which Hagedorn predicts to be independent from the electron-positron collision energy W in 
the mass center referential, and c is a constant. This theory fits the data quite well for small 
W, say W < 10 Gev, but exhibits a pronounced failure for W increasing up to say 160 Gev. 
Very recently, Bediaga, Curado and Miranda [|102|| have used, along Hagedorn's lines, the 



present generalized statistics. The results are indicated in Figs. 8 and 9. Remark that (i) 
q varies smoothly and monotonically with varying W (Hagedorn's theory is recovered in 
the W — ► limit), and (ii) To ~ 0.11 Gev and practically independs from W as desirable 
from Hagedorn's arguments. These results can be considered as a strong evidence of the 
applicability of the nonextensive thermostatistics to specific anomalous systems. 

G. Emulsion chamber observation of cosmic rays 

Cosmic rays can be observed by using a variety of detectors (such as Pb detectors; see 
|103|| and references therein). Typically, showers of (clustered or individual) elementary 



particles appear which start at the so called vertex. These vertex are localized at various 
depths. The distribution of their depths can be measured (see and references therein for 



the measurements done at the Mount Pamir lead chambers) . This distribution was recently 
fitted by Wilk and Wlodarcsyc | |105| | with the q — 1.3 function which emerges within the 
present formalism. 



H. Reassociation of heme-ligands in folded proteins 

In the folded conformational state, proteins might exhibit fractal effects. One such case 
might be the time evolution of the re-association of molecules that have been taken away 
from their natural positions. For instance, if O2 molecules are dissociated, through light 
flashes, from their natural Fe positions in a heme protein and reach positions outside the 
heme pocket, they tend to start rebinding, and, for so doing, they might have to follow 
a fractal path, or be under the dynamical influence of fractal excitations (e.g., fractons). 
Anyhow, this re-association phenomenon has been lengthily studied by Frauenfelder et al 
||106|| . If we define £ = N(t)/N(0) where N(t) is the number of molecules that have not 
yet re-associated at time t, the monotonically vanishes with t. The results obtained by 
photo-dissociating CO molecules from Sigma Type 2 sperm whale Myoglobin (Mb) dissolved 
in a glycerol- water solution are shown in Fig. 10. For times not too long, the experimental 
data have been fitted by Frauenfelder et al ||106|| with 

i = {l+t/to)- n (88) 

where to and n smoothly depend on the temperature T. Bemski, Mendes and myself |107|| 
have argued that, within the generalized formalism, the following equation naturally appears: 

§ = -A, £ 9 {X q > 0; q > 1) (89) 
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Its solution is given by 

1 



[l + (q-l)X q t] — 



(90) 



This expression recovers, for q — 1, the usual exponential relaxation, and reproduces the 
Frauenfelder form through the identifications l/(q — 1) = n and l/[(q — l)A g ] = t . Besides 
reobtaining the Frauenfelder empiric law, the present scheme allows for a better approxi- 
mation if a crossover is admitted. More precisely, the above differential equation can be 
generalized as follows: 

§ = -Mrf-(^-^H (r<q) (91) 



The general solution involves ||107|| a hypergeometric function. The fitting is shown in Figs. 
10 and 11. A detailed model which would justify the above phenomenological differential 
equation would be welcome. 



I. Diffusion of Hydra Vulgaris 

Upadhyaya et al ||108|| are presently performing interesting experiments on Hydra Vulgaris 
(a cylindrical body column with inner and outer cells, respectively referred to as endodermal 
and ectodermal respectively) in physiological solution. The endodermal cells are more adhe- 
sive than the ectodermal ones. The authors have measured the velocity distribution PdV^I) 
of the "vertical" component of the velocity during the diffusion of endodermal Hydra cells 
in an ectodermal aggregate. The results are presented in Fig. 12, where the velocity unit is 
10~ 6 m/hour and the probability is represented by the histogram of the number of counts. 
These results were fitted with 

"<W> - (TTW (92) 

with the values of (a, b, c) indicated in the figure. Through the identification 

a = P(0); b=(q-l)/V 2 ; c = (93) 

we precisely have the law which emerges within the present formalism, namely 

Pm) = [i + (g-i)(VW/ (<? - 1} (94) 

with q = 1.53. The next desirable step of course is to formulate a specific model for Hydra 
which would lead to this law, but this remains to be done. 
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J. Citation of scientific papers 



An interesting study was recently done by Redner ||109j ], in which the statistics of citations 
of scientific papers is focused. He exhibited the number N(x) of papers which have been 
cited x times for two long series, namely one (6 716 198 citations of 783 339 papers) from the 
Institute of Scientific Information (ISI) and another one (351 872 citations of 24 296 papers) 
from the Physical Review D (PRD). As expected, in both examples, N(x) monotonically 
decreases with x. Redner fitted the (relatively) low-x data with a stretched exponential of 
the form 

N(x) = N(0) e - {x/xo)0 (95) 

with (3 = 0.44 and 0.39 for the series ISI and PRD respectively. Also, he remarked that the 
large- x data exhibit a power law, namely close to oc 1/x 3 . He argues that this different 
functional behavior for low and large values of x must reflect different phenomenologies in 
these two regimes. In contrast with this viewpoint, Albuquerque and myself ||110|| argue 



that this is not necessarily so since the data can be quite satisfactorily fitted with a single 
function, namely 

M x ) = t , N ^\] , n n (96) 

V ; [1 + (q - l)Xx]9/(fi-i) v ; 

with q = 1.53 and 1.64 for the series ISI and PRD respectively: see Fig. 13 The satisfactory 
quality of the fittings is, after all, not so surprising, since we have mentioned earlier in this 
paper the connection [S4j] of this formalism with the Zipf law. 

K. Electroencephalographic signals of epilepsy 

It is since long known that the analysis of signals can be done within formalisms which 
use entropic forms. One such application has been recently done on EEG records of epilleptic 



humans and turtles ||111|| . The simultaneous use of wavelet-based multiresolution analysis 
including the nonextensive entropy S q leads to signals whose interpretation can be clini- 
cally neat and pharmacologically convenient. The authors of this novel processing suggest 
perspectives for building up automatic detection devices. 

L. Cognitive psychology 

The development of artificial neural networks and their connections with statistical me- 
chanics (e.g., the Hopfield model for associative memory) makes quite natural the approach 
of cognitive problems with the present nonextensive formalism. Within this philosophy, we 
performed ||112|| an experiment of learning/ memorization (of 5 x 5 and 7x7 square matrix 



having circles and crosses randomly distributed once for ever) with students of the Univer- 
sity level; 150 students were interviewed, the first 30 in order to optimize the experimental 
protocole, and the other 120 to make the measurements of the time-evolution of the total 
amounts of errors when the original matrix was successively shown and hidden. The average 
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results were then fitted with those obtained, for the same task, with a learning machine ||113| 



having a perceptron architecture and an internal dynamics based on the Langevin equation 
|3| generalized by Stariolo to arbitrary q. The (average) learning time of the machine 



turned out to monotonically increase with q, exhibiting a practically divergent derivative 
at q — 1. The best human- machine fit occurred for q slightly above unity More experi- 
ments and comparisons along these lines would be very welcome. Indeed, they would help 
better understanding some cognitive phenomena, on one hand, and could generate efficient 
machines for specific tasks, on the other. 

M. Turbulence and time evolution of financial data 



In 1996 Ghashghaie et al ||114|| compared financial data with those obtained from turbu- 



lent behavior and showed very similar behaviors when appropriate scalings are used. Ramos 



et al ||115|| have recently shown that all these data can be satisfactorily fitted with the 
functional form which emerges from the present formalism. Olsen and Associates data con- 
taining bid-ask quotes for US dollar-German mark exchange rates (1,472,241 records) are 
presented in Fig. 14 (probability density P^t^An) of price changes; An = n(t) — 7r(t + At) 
with At = 640s, 5120s, 40960s, 163840s from top to bottom in the figure). The turbulent 
flow data [ |116|| are presented in Fig. 15 (probability density P& r (Av) of velocity differences; 
Av = v(r) — v(r + Ar) for spatial scale delays Ar = 3.3r], 18.5r], 138^, 325r] from top to 
bottom in the figure, where rj is the Kolmogorov scale, i.e., the critical limit for occurence 
of viscous dissipation). All these data exhibit a slight left-right assymmetry, which has been 
taken into account by Ramos et al: they used the same q for both sides but different widths. 
Needless to say that specific models leading to these fitting functions are very welcome. 

V. COMPUTATIONAL EVIDENCES AND CONNECTIONS 
A. Thermalization of a hot gas penetrating in a cold gas 



In 1991, Waldeer and Urbassek ||117|| made, assuming d = 3 Newtonian mechanics, a 
computational simulation in which a certain amount of high energy particles penetrate into 
a cold gas and are thermalized through the interactions between molecules. The cold gas is 
initially put at BG thermal equilibrium at temperature Tc- The high energy particles at time 
t = are randomly distributed in energy at a quite high energy per particle. The interaction 
potential was assumed to be hard sphere at short distances and decreasing, at long distances, 
like r~ a . They analyzed three typical situations, one with a — * oo, hence well above d (i.e., 
very short range interactions), the second one with a = 4 (i.e., short range interactions), 
and the last one with a = 8/3, which is below d (i.e., long range interactions). In their 
simulation, they follow the time evolution of the energy distribution of the hot particles. 
After a transient, this distribution evolves with a regular pattern. For a > d, this pattern 
basically is the BG distribution with a temperature T(t) which gradually approaches Tc 
from above (with Um t ^ooT(t) = Tc), in other words, through curves which approximatively 
are straight lines in a log-linear plot. For a < d, this approximation occurs through curves 
which are close to straight lines... in a log-log plot ! (Notice that the curvature in log- log 
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plots tends to be upwards for a < d, whereas it is downwards for a > d; see Figs. 1, 2 and 
3 of ||117| ). This power-law behavior is typical of q > 1. This pecualiarity was invoked by 



Koponen || in 1997 as a justification for using the present generalized formalism to discuss 
electron-phonon relaxation in ion-bombarded solids if the interactions are long-ranged. A 
study like that of Waldeer and Urbassek [ |117]| which would systematically address the details 
of that thermalization by gradually varying a across d is missing and would certainly be 
very welcome. 



B. Long-range classical Hamiltonian systems: Static properties 

Let us focus here on what we refer to as weak violation of BG statistics. We use this 
expression to distinguish it from what we call strong violation of BG statistics. Both of 
them lead to nonextensive quantities, but, whereas the strong violation concerns q ^ 1, 
the weak one concerns q = 1 calculations. To make all this explicit we shall here focus on 
classical systems, i.e., all observables are assumed to commute. Let us consider the following 
paradigmatic Hamiltonian: 

# = if>*+E^) (97) 

i=i i^j 

where m is a microscopic mass, are the (i-dimensional linear momenta and positions 

associated with N particles, and = rj — r». A typical situation is that of a finite confined 
system but, if some care is taken, the system could as well be thought of as having periodic 
boundary conditions. To be specific, let us assume 

V{r l3 ) = 4l-^ (A>0,5>0,0<a<12) (98) 

r ij r ij 

where, in order to avoid any singularity at the origin (for any dimension d not exceedingly 
high), we have assumed, for the repulsive term, the Lennard- Jones exponent 12. What we 
desire to focus on in the present discussion is possible singularities associated with infinite 
distances, i.e., the effects of long- range (attractive) interactions. The case (a, d) = (6,3) 
precisely recovers the standard Lennard- Jones fluid; the case (a, d) = (1,3) is asymptot- 
ically equivalent to Newtonian gravitation; the case (a, d) — (d — 2, d) is asymptotically 
equivalent to (i-dimensional gravitation (i.e., the one associated with the solutions of the 
rf-dimensional Poisson equation); the case (a, d) = (3,3) basically reproduces the distance 
dependance of permanent dipole-dipole interaction. The range of the (attractive) interaction 
increases when a decreases; a — > 12 corresponds to very short-ranged interactions, whereas 
a = corresponds to the situation of the Mean Field Approximation, where every particle 
(attractively) interacts with every other with the same strength, in all occasions. 

A typical quantity to be calculated within BG statistics is the following one (basically 
related to the T = internal energy per particle): 



oo 

d-1 



dr r d ~ l r~ a (99) 
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where the distances r have been expressed in units of a characteristic length of the problem. 
We immediately verify that this integral converges if a > d, and diverges if < a < d. 
Consequently, thermodynamic calculations in the < a < d case have to be done with some 
care, and not blindly following the standard rules associated with BG statistics (i.e., q = 1). 
It is in this sense that we use the expression "weak" violation of BG statistics. The care 
to which we refer is the fact that we have to strictly consider the finite size of the physical 
system. Consistently, a relevant quantity that emerges naturally is 

N* = d dr r d ~ l r- a = — (100) 

Ji 1 — a/d 

We can check that, in the iV — > oo limit, we have 



i 



if a/d > 1; 



a/d-l 

N* ~ { In N if a/d = 1; (101) 
if 0<a/d<l. 



1-a/d 



As it will become transparent later on, what these regimes imply is that the system is 
extensive for a/d > 1 (hence standard thermodynamics apply), whereas it is nonextensive 
for < a/d < 1, and special scalings become necessary [ |118|jll9|| in order to have both 
a mathematically well posed problem, and a physical unfolding (or qualification) of the 
nonextensive region. The a/d > 1 regime has since long been analyzed |120|| , and it is 
well known that extensivity (or, stability, as also referred to) is lost for a/d < 1. However, 
to the best of our knowledge, the scalings associated with the quantity iV*, as well as its 
numerically efficient collapsing properties, were introduced for the first time by Jund et al 
[Tn§ in 1995. 

A quantity related to N*, namely N turns out to be even more convenient. It is defined 



through [121 



iN^^-a/d] , s 

N = N + 1 = - — ^ 102 

1 — a/d 



In the N —>■ oo limit, we have 



if a/d > 1; 

N ~ { In N if a/d = 1; (103) 
if 0<a/d<l. 



1-a/d 



In the limit a/d —>■ oo, N —>■ 1; in the limit a/d — > 1 + 0, iV~ l/(a/d — 1); in the limit 
a/d — > 0, iV ~ N. Roughly speaking, iV characterizes the effective number of neighbors 
that can be associated with a given particle. This is the convenience to which we referred 
above. 

We are ready now to present the kind of size-scalings we expect to be necessary for 
thermodynamically describing a generic classical Hamiltonian system with the type of inter- 
actions above mentioned. Let us focus on a simple fluid, and start with the standard case, 
i.e., a > d. Its Gibbs energy G(T,p,N) is given by 
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G(T, p, N) U(T, p, N) S(T, p, N) V(T, p, N) 

N N N N 1 1 

where U, S, V, N, T and p respectively are the total internal energy, total entropy, total 
volume, total number of particles, temperature and pressure. In the N — > oo limit, we 
obtain 

g(T,p) = u(T,p) - T s(T,p)+pv(T,p) (105) 

where the corresponding densitary variables have been introduced. 

In contrast with the above, if we have < a < d, the scalings are different, namely 

G(T,p,N) ^ U(T,p,N) _ T S(T,p,N) p_ V(T,p,N) 
NN ~ NN N N N N 

Consistently we have 

g(f,p) = u(f,p) - f s(T,p)+pv(T,p) (107) 

where 

f=t; p=l (108) 
N N 

These equations recover the previous ones, i.e., those associated with the a > d, as a 
particular case. Indeed, for a > d, N becomes a constant. 

So, we see that long range interactions have important thermodynamical consequences, 
namely 

(i) the energy quantities (G, U, which normally appear alone) that were extensive for 
a > d loose their extensivity; 

(ii) the non-energy quantities (S, V, which normally appear in canonical pairs with 
intensive quantities) that were extensive for a > d preserve their extensivity; 

(iii) the control parameters (T, p) that were intensive for a > d loose their intensivity. 
Consistently, to have mathematically well defined and physically useful equations of 

states and related quantities, everything must refer to finite quantities, hence, we must 
express all relations with the above rescaled variables. This does not imply that thermal 
equilibrium occurs through sharing equal values of T, p, etc. The zero-th principle of 
thermodynamics appears to hold in the usual way, even if we have long range interactions 
in the system. Although we have illustrated these features on a fluid, it is clear that the 
same considerations hold for all types of thermodynamical systems (magnets, dielectric 
substances, elastic solids, etc). 

Two particular remarks must be made at this point: 

(i) When every element of the system equally interacts with each other (i.e., a — 0), 
N = N, and consequently T — T/N. In what concerns the thermostatistical approach of 
a system, this is equivalent to dividing the microscopic coupling constants by N, a familiar 
feature that is artificially imposed in all Mean Field calculations. We have used the word 
"artificial" because, whenever a < d, the Hamiltonian which includes the microsopic inter- 
actions indeed is nonextensive! (and so is U). To divide the two-body coupling constants by 
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N when a = (or, by N, when a > 0) certainly is an artificial manner of forcing to be ex- 
tensive a Hamiltonian which physically is not. This practice is traditionally frequent among 
magneticiens (who divide J by N), but certainly not among astronomers, who normally do 
not even think about what would be a very strange way of renormalizing the gravitational 
constant G ! 

(ii) If a singularity (for example, a critical phenomenon) occurs under particular physical 
conditions, it must generically occur at finite values of T, p, etc, and not at finite values 
of T, p, etc. Let us illustrate this for the simple case of a critical temperature: T c must 
be finite, hence, if a/d > 1, it must be T c (a,d) oc l/(a/d — 1) This implies that T c must 
generically diverge for all classical systems at the extensive-nonextensive frontier. More 
precisely 

T c (a,d) ^p- (a/d ^1 + 0) (109) 

a/d — 1 

where A(d) is a system-dependent finite constant. In fact, let us anticipate that this precise 
behavior has been observed in the systems available in the literature [|119|| (d = 2 and d = 3 
Lennard- Jones-like fluids, d = 1 and d = 2 Ising and Potts ferro- and antiferromagnets, 
etc) , with no exception. To illustrate the connection between fluid models like the extended 
Lennard- Jones one above considered, and localized spin systems, let us briefly focus on the 
Ising ferromagnet. The simplest long-range iV-spin cubic-lattice Hamiltonian of this kind is 
given by |T22] 

h = ( J > °; a ^ °; s i = ±x yi ) ( n °) 

where, for d — 1, =1, 2, 3, for d = 2, rij = 1, 2 1 / 2 , 2, for d = 3, 



r 



1, 2 1 / 2 , 2 1 / 3 , 2,... ; and so on. Clearly, the limit a — >• oo recovers the first- neighbor 
c?-dimensional spin 1/2 ferromagnet, whereas a = corresponds to the Mean Field Approx- 
imation. For this model, kTc(a, d) / J diverges for < a/d < 1, decreases for a/d increasing 
above unity, and approaches the first-neighbor value (0 for d — 1, 2.269... for d = 2, etc) for 
a/d — > oo. Also, kT c / J ~ A((i)/[(a/(i) — 1] in the a/d ^1 + limit. The introduction of 
T nicely enables the unfolding of the region where T c diverges. Indeed, kT c /J = kT c /(JN) 
is finite in both extensive (a/d > 1) and nonextensive (0 < a/d < 1) regions, thus providing 
an enlightening unification. Finally, let us mention that, it seems that all equations of states 
(e.g., limN^ 00 M(T 1 N)/N can be, for all a/d < 1, mapped into that associated with the 
Molecular Field Approximation. This simplifying feature appears to hold only for the static 
thermodynamic properties, and not for the dynamical ones, as will become clear later on (in 
Subsection V.G). 

To close this subsection, let us emphasize that what we have been focusing on here is 
what we refer to as the weak violation of BG statistical mechanics (see Fig. 4). These 
are analytical or Monte Carlo q = 1 calculations (i.e., the energy distribution obeys the 
Boltzmann factor), but the variables must be scaled with N, which is not at all necessary 
for the standard, short-ranged interacting systems. 
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C. Long-range tight-binding systems 



In the previous subsection we addressed classical systems. It is clear, however, that 
similar nonextensivity is expected to emerge in quantum systems if long-range interactions 
are present. One such Hamiltonian is the tight-binding-like which follows [ 123, 124 1 : 



H 



N 
i=l 



V 



111) 



where cf and q are the creation-annihilation operators associated with electrons on site i, 
the {ej} are the on-site energies, and V is the inter-site energy. The T = electron diffusive 
properties corresponding to this Hamiltonian exhibit a variety of anomalies intimately re- 
lated to N, as preliminary shown by Nazareno and Brito | ]123}| and studied with more details 
in PI. 



D. Granular systems 

In 1995, Taguchi and Takayasu || simulated a vertically vibrated bed of powder with 
inelastic collisions and studied the distribution of horizontal velocities. In the lower layers (so 
called solid phase) they observed a standard Maxwellian (Gaussian) distribution of velocities. 
The situation was sensibly different in the upper layers (so called fluidized phase). Indeed, 
there the distribution was a Student's t-distribution, precisely the one appearing in Eq. 
(42) with q = 3, and assuming an energy proportional to the (velocity) 2 (together with 
d = r = 2, hence a constant density of states). This anomaly must be related to the fractal- 
like granular clusterization which occurs in real space [|125|| but a deep analysis would be 



welcome. Also would further simulations, for instance of the cooling type. Studies of such 
computational models, either externally forced or just left to their own isolated evolution, can 
provide important physical insights, especially if quantities like the energy distribution, the 
Lyapunov spectrum (or at least its maximum value) or possible multifractality are focused 
on. 



E. d = 1 dissipative systems 

One-dimensional maps constitute the simplest systems which might present chaos. Ba- 
sically they consist of the following recurrent equation: 

x t+1 = h(x t ; a) {t = 0,1,2,...) (112) 

where h(x; a) is a rather simple nonlinear function, and a is a control parameter. Typically, 
both x and a are real numbers (but higher- dimensional situations are of course possible, 
and frequently studied). The logistic map, for instance, exhibits this structure. Typically, 
for a < a c , the system exhibits simple orbits, the attractor being a cycle whose number of 
elements is finite. For a > a c , the system can exhibit attractors with an infinite number 
of elements. The value a c is the critical one, usually referred to as the chaos threshold; 
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the associated attractor typically constitutes a multifractal characterized by the so called 
multifractal function f(an), where ajj is the Holder exponent. The function is 

generically concave, attains its maximum at a value of an in the interval [a™ n , a™ x \ and 
this maximum equals the fractal or Haussdorf dimension df. An important feature of this 
type of maps is the sensitivity to the initial conditions (and, of course, the rounding at any 
intermediate calculation). More precisely, if we note Ax(0) a small variation in the initial 
condition xq, and follow its time evolution Ax(t), we can define the sensitivity function 
as £(t) = ^mAa;(o)->o Ax(o) ' ^ mos ^ values of a, satisfies d^/dt = Ai£, hence 

£ = expAi£ (113) 

where Ai is the so called Lyapunov exponent. If Ai > (Ai < 0) the system is said strongly 
sensitive (insensitive) to the initial conditions. The Ai = possibility can also occur and is 
referred to as the marginal case. In this situation f (t) satisfies PL[126i[l27| d^/dt = X q £*, 
hence 

£ = [l + (l-g)A,t]^? (114) 

which recovers Eq. (25) as the q — 1 case. Two Ai = possibilities exist, namely q < 1 with 
X q > (weakly sensitive to the initial conditions), and q > 1 with X q < (weakly insensitive 
to the initial conditions). For instance, the logistic map exhibits q — 1 for almost all values 
of a but exhibits q < 1 at the chaos threshold and q > 1 at every doubling-period as well 
as tangent bifurcations. Moreover, it has been shown that a large class of such systems (for 
which f(ajj in ) = f(<y 7 fi ax ) = 0) verify, at the chaos threshold, the following scaling law ||127|| : 

i i i 

— - — (115) 



1 - g ag* 



This is a fascinating relation. Indeed, its left-hand member concerns the dynamics of the 
sensitivity to initial conditions of the map, whereas its right-hand member concerns pure, 
though nontrivial, geometry. Under what precise mathematical conditions does it hold? 
How should it be generalized in order to also cover the standard case of Euclidean ge- 
ometry (a™ n = a^ ax = df = 1) for which one expects q = 1? (Should we also con- 
sider simultaneously /(a™ n ) = f(ot max ) = 1 ?; Could Eq. (115) be generalized into say 
1/(1 — q) = l/[«™ n — /(a™ n )] — l/[a]} ax — f(a I } ax )}, the Euclidean case thus corresponding 
to a special limit of the type ot™ n / f(a™ n ) = a I r } ax /f(a]} ax ) = q = 11). What happens for 
two- or more-dimensional maps? What happens if, instead of maps, we have ordinary (or 
even partial) differential equations? To answer all these questions, computational effort is 
invaluable. 

Before closing this Subsection it is mandatory to clarify what the index q appearing in 
the differential equation yielding Eq. (114) has to do with the one appearing in the present 
generalized entropy. In fact, they are one and the same, and the connection is established 
through the so called Pesin equality or identity. Let us illustrate the basic ideas on the 
logistic map herein considered. Assume that we make a partition of the x interval into 
a large number M of equally small windows, chose arbitrarily one of those windows and 
randomly put a large number iV of points inside. We then calculate (by using the set of M 
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probabilities corresponding to the ratios of numbers of points belonging to each window) the 
t = value of the BG entropy Si(0), which is going to be very close to zero (strictly zero in 
the (M, N) — > (oo, oo) limit). We then allow each of the iV points to evolve according to the 
logistic map until an attractor is achieved. The entropy Si(t) will grow with t until arrival 
to a saturation value ^(oo) which depends on (M, N) (necessarily Urn tv_ >00 S'i(oo) < In M). 
In the M — > oo limit, the growth of S\ (t) is in fact linear (see, for instance, |[L28|| , which 
enables the following characterization of the so called Kolmogorov- Sinai entropy: 

Si(t) 

Ki = limt^oolimM^oolirriN^oc-^ — (116) 

Quite generically, the Pesin inequality holds, which states (for one-dimensional nonlinear 
dynamical systems) that 

Ki < Ai if Xi > (117) 
and K\ = otherwise. Since we are only presenting a sketch of the situation, let us from 



now on address those particular systems for which the equality holds ||129|| . For those it is 

K x = Ai if Ai > (118) 

and Ki = otherwise. This type of analysis is convenient either if we have simple orbits (i.e., 
strong insensitivity to the initial conditions, i.e., for Ai < 0) or if we have strong chaos (i.e., 
strong sensitivity to the initial conditions, i.e., Ai > 0). But this analysis is a very poor one 
at say the edge of chaos, where Ai = and we have weak sensitivity to the initial conditions. 
It is to unfold this type of situation that S q becomes extremely useful. Let us show how. 
At the chaos threshold we have K\ = Ai = 0. But if we follow the same procedure we just 
described for calculating K\, but using instead S q (t), an interesting phenomenon can be 
revealed, which we describe now. Let us first define the following generalized Kolmogorov- 
Sinai entropy: 

S (t) 

K q = lim t ^ c Jim M ^ c Jim N ^- 1 Y- (119) 



A value q* is generically expected to exist ||130|| such that, for q > q* (q < q*), K q = 



[K q diverges), and K q * is finite ! Furthermore, it can be argued [Q that the above Pesin 
equality can be generalized as follows: 

K q = \ if\ q >0 (120) 

and K q = otherwise. It is through this important type of (in)equality that the connection 
emerges between S q and the power-law time-dependence of the sensitivity to the initial 
conditions. The particular value q* above described is what was numerically calculated in 
38| , | 1 2 6| , p~2T|| , and satisfies the scaling (115) 



Some of the above statements can be trivially checked with the logistic map at its chaotic 
region (i.e., for Ai > 0). We know in that case that Si(t) oc t, hence (assuming the simple 
case of equiprobability) the total number of possibilities W(t) grows exponentially with t. 
For any q > 1, S q (t) is always bounded, then K q necessarily vanishes. For any q < 1, S q (t) 
grows like the 1/(1 — q) power of W(t), which in turn, as said before, grows exponentially 
with t, hence necessarily K q — > oo. We conclude that q* = 1 ! The same picture is expected 
to hold for weak chaos. 
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F. Self-organized criticality 

In the previous example, fine tuning (e.g., a = a c ) is necessary to observe the anomalous 
(q 7^ 1) behavior. Let us address dissipative systems with many degrees of freedom, very 
particularly those which do not need fine tuning. Would robust systems like those exhibiting 
self-organized criticality [|131|| (SOC) also present q 7^ 1 behavior? The answer is yes, as it 
has been clearly exhibited in at least three computational systems, namely the Bak-Sneppen 
model for biological evolution, the Suzuki-Kaneko model for imitation games and the Bak- 
Tang-Wiesenfeld model for sandpiles [|132|| . In these systems, the Hamming distance plays 
the role played by £ in the previous ones. Also, the relevance of the order of the t — > 00 
and N 00 limits has been exhibited. Like in the conjectural Fig. 4, the q 7^ 1 behavior is 
observed only in the limt^ool^N ~^oo order. On what model ingredients does q depend? Is a 
taxonomy in universality classes analogous to that of standard critical phenomena possible? 
Is a multifractal f{a.L) function hidden somewhere? Does a scaling law like that of Eq. (115) 
still hold? Again, additional computational effort is very welcome. 



G. Long-range classical Hamiltonian systems: Dynamic properties 

Let us finally focus on the "heart" of statistical mechanics, the dynamics of the systems 
on which Boltzmann himself was meditating, namely the Hamiltonian systems with many 
degrees of freedom. Although lots of interesting quantum nonextensive phenomena must 
exist, here we shall restrict ourselves to the classical ones. We expect them to be able 
to provide nonextensive anomalies in a kind of pure, or simpler manner. Since a classical 
canonical Hamiltonian must satisfy the Liouville theorem, the Lyapunov spectrum must 
be symmetric with respect to zero, the corresponding eigenvalues being necessarily coupled 
in pairs of positive and negative values with the same absolute value. Consequently, to 
study the sensitivity to the initial conditions it suffices to study the maximum Lyapunov 
exponent. If it is positive, the system will generically be strongly chaotic, and will therefore 
easily satisfy the ergo die /mixing hypothesis (equality of time and ensemble averages). If, 
however, the maximum Lyapunov exponent vanishes, the entire spectrum will necessarily 
vanish, hence the system will be, at most, weakly chaotic, and will therefore have difficulties 
in satisfying the ergo die /mixing hypothesis, at least at not extremely large times (reflecting 
the macroscopic size of the system). The d = 1 coupled planar rotators iV-body model with a 
two-body coupling constant proportional to l/r a (r = distance between two given rotators) 
has been recently studied (for a = in ||133 |, and, for a > in ||134|| ) in the microcanonical 



ensemble. It has been established that, above a critical (conveniently normalized) total 
energy, the maximum (conveniently normalized) Lyapunov exponent is, in the N — ► 00 
limit, positive (zero) for a > 1 (0 < a < 1). More precisely, this maximum Lyapunov 
exponent is proportional to 1/N K where k{q) appears to be a monotonic function which 
decreases from k(0) to zero while a increases from to 1, and remains zero for all a > 1). It 
must be recalled that it is only for a > 1 that the standard BG prescriptions provide finite 
integrals in the relevant calculations. If we were to discuss the d-dimensional version of the 
same model, we would certainly have n(a, d). It is certainly possible that it is n(a/d), and it 
would be interesting to check such a hypothesis. If we were to consider not planar (like the 
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XY ferromagnet) but rather three-dimensional (like the Heisenberg ferromagnet) rotators, 
would k be insensitive to that, or it would depend on the specific model? What would 
happen if, instead of localized rotators, we were to consider a long-ranged Lenard- Jones-like 
fluid, or d- dimensional gravitation? All these questions are certainly interesting, and worthy 
of further computational efforts. 



In a recent paper, an essentially a = model was considered ||135|| , and, under certain 
circumstances, a crossover was found between anomalous (at times smaller than t(N)) and 
normal (at times larger than t(N)) diffusion, with t(N) oc N. What happens if a > 0? 
Does r scale like N = N*(l + a/d) = (N l - a ' d - 1)(1 + a/d) /{I - a/d) ? What happens for 
other models? The behavior observed for the particular model that was studied is consistent 
with Fig. 4. But is it exactly that conjecture that is going on? Only the study of the energy 
distributions (of single particles or of relatively large subsystems of the iV-body system) 
themselves can provide the answer. What about the distributions of velocities? Are they 
Maxwellian (i.e., Gaussian) for a/d > 1 and non-Maxwellian otherwise? Are they Levy's 
or Student's t-distributions for a/d < 1? If so, what is the dependence of q{a,d) ? Maybe 
q(a/d) ? Are the associated fluctuations anomalously time-correlated? Can nonmarkovian 
processes be present when the system is nonextensive (i.e., when < a/d < 1)? Plenty 
of intriguing questions that, sooner or later, will have to be answered, mainly through 
computational work (at least the first approaches). Better sooner than later! 



H. Optimization techniques; Simulated annealing 

The so called Optimization problem consists basically in determining the global minimum 
(or minima, if degeneracy is present) of a given cost funtion E(x), where x is a discrete or 
continuous (^-dimensional variable. This problem can become extremely complex depend- 
ing on the cost function having a large number of local minima, and on the dimension d 
being high. For the ubiquitous cases (in physics, chemistry, neural networks, engineering, 
finances, etc) for which analytic discussion is not tractable, a variety of computational al- 
gorithms have been developed. A special place among these is occupied, because of its 
efficiency and paradigmatic value, by the Simulated Annealing (SA) introduced in 1983 by 
Kirkpatrick et al [|136|| . Its denomination comes from its total analogy with the well known 
annealing technique, frequently used in Metallurgy for making a molten metal to reach 
its crystalline state (global minimum of the relevant thermodynamic energy). In SA, one 
or more artificial temperatures are introduced and gradually cooled, acting as a source of 
stochasticity, extremely convenient for eventually detrapping from local minima. Near the 
end of the process, the system hopefully is in the attractive basin of one of the global minima, 
the temperature is practically zero, and the algorithm asymptotically becomes a steepest 
descent one. The challenge is to cool the temperature the quickest we can but still having 
the guarantee that no definitive trapping in any local minimum will occur. More precisely 
speaking, we search for the quickest annealing (i.e., in some sense approaching a quenching) 
which preserves the probability of ending in a global minimum being equal to one. SA 
strictly follows a BG scheme. Let us illustrate for continuous x. The system "tries" to visit, 
according to a visiting distribution assumed to be Gaussian in the neighborhood of its actual 
state. The jump is always accepted if it is "downhill", i.e., if the cost function decreases. If 



35 



it is "uphill", the jump might be accepted with a probability given by the Boltzmann factor 
corresponding to that cost function. It has been shown that the probability of ending on the 
global minimum equals unity if T(t) decreases logarithmically with time t. This algorithm 
is sometimes referred to as Classical Simulated Annealing (CSA) or Boltzmann machine. 
We easily recognize that, if instead of decreasing, the temperature was maintained fixed, 
this procedure precisely would be the well known Metropolis et al one for simulating BG 
thermostatistical equilibrium. 

This optimization machine has been generalized within the present statistics as follows 
|58fl . The visiting ditribution is generalized to be a qy-Gaussian, and the acceptance Boltz- 
mann factor is generalized to be a g^-generalized factor, where qy and respectively are 
the visiting and acceptance entropic index. The cooling schedule is generalized as follows: 

This is the Generalized Simulated Annealing. This machine is characterized by (qv,QA)- 
The choice (1, 1) corresponds to CSA, and the choice (2, 1) corresponds to the so called Fast 
Simulated Annealing (FSA). The CSA corresponds to a cooling given by 

T(t) = T(l) — ^ (t = l, 2, 3, ...) (122) 

The FSA corresponds to a faster cooling given by 

T(t)=T(l)i (t = l, 2, 3, ...) (123) 

The limiting case qy = 3 corresponds to 

T W = T W (t + l 3 )2-l (t = X ' 2 ' 3 ' - ) (124) 

These particular cases illustrate the great computational advantage that can be obtained 



by speeding up the algorithm by conveniently choosing qA (see also |68fl). In practice, a 



convenient choice for qA is slightly below 3. The choice of qy seems to be more model- 



dependent. Details can be seen in a by now vast literature |p9|-pq], in which applications 
have been done and variations have been performed concerning a variety of classical and 
quantum physical problems, the Traveling Salesman Problem, and many others. The first 
application |)7]] in quantum chemistry concerned simple molecules of the series CH 3 — R 
and some others, including the H2O3 one, by using the MOPAC program package. Nowa- 
days, Straub (in Boston), Okamoto (in Okazaki), and Ellis-Mundim-Bisch (respectively in 
Chicago, Salvador and Rio de Janeiro) and co-workers are currently improving and applying 
these techniques to complex molecules such as polypeptides, in particular with the aim of 
studying the important, though hard, protein folding problem. 



VI. FINAL REMARKS 

Boltzmann-Gibbs statistical mechanics and standard thermodynamics do not seem to be 
universal. They have domains of applicability quite poorly known nowadays. The precise 
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knowledge of the restrictions for their validity is conceptually and practically very impor- 
tant. A nonextensive generalization of these formalisms is now available PiPl; see Table I. 
It has been developed to cover at least some of these difficulties. Several important types of 
systems have been focused on in the present paper, which should substantially clarify the 
situation. These efforts span a wide epistemological variety, which goes from clear-cut theo- 
ries to phenomenological ones, to quite well or less well understood fittings and connections. 
As exhibited at length here, the areas on which this formalism has been satisfactorily ap- 
plied includes physics, astronomy, chemistry, mathematics, biology, economics, linguistics, 
cognitive psychology, etc. However, in spite of all this sensible progress, some inter-related 
crucial points are still to be understood and established on a neat and transparent basis. 
These include (i) the zeroth principle of thermodynamics and its connections with the ther- 
modynamic limit, properties which would in principle exhibit the mathematical connection 
of the weak violation of the BG statistics (i.e., introduction of N within the q = 1 formal- 
ism) with the strong one (q ^ 1); (ii) the functional dependence of q on (a, d) for long-range 
interacting Hamiltonians (0<a/d<l) and its connection with anomalous diffusion; (iii) 
the physical interpretation of the {pi} distribution and of the escort one {Pi}, as well as 
clear-cut prescriptions for using one or the other when fitting experimental data (in the 
meanwhile, it appears that {Pi} is the one to be used for equilibrium distributions of ki- 
netic, potential or total energies, whereas {p^ is the one to be used for real-space diffusive 
nonequlibrium phenomena, either Levy or Student's [|137|| distributions according to whether 
correlations between jumps are suspected to be absent or present); (iv) what are the generic 
physical conditions for using the microcanonical, canonical and grand-canonical ensembles, 
under what exact and fully specified conditions they are expected to be thermodynamically 
equivalent, and the possible relevance for the so called thermogravitational catastrophe; (v) 
the clear connection with microscopic dynamic properties such as (partial) lack of ergodicity 
and mixing, the generalization (to weak chaos) of the Pesin inequality, the complete domain 
of validity of scaling relations connecting q to multifractality, and the possible relevance for 
SOC, spin-glasses and similar phenomena; (vi) the clear physical connection with quantum 
groups and, in general, deformations of relevant Lie algebras, and through these, the possible 
relevance for quantum gravity and the deep (possibly discrete, multifractal- like) structure 
of space-time. 

On speculative grounds, one might think of two conjectures, to be clarified (i.e., rigorously 
formulated), confirmed or refuted. The first of these conjectures can somehow (on intuitive 
grounds) be formulated as follows. Strongly ergo die /mixing phenomena are ubiquitous in 
Nature; essentially, they are driven by microscopic interactions which are short-ranged in 
space-time (short-range forces, short-range memory, nonfractal boundary conditions); their 
basic geometry tends to be continuous, Euclidean-like; their thermodynamics is extensive; 
their central laws (energy distribution at equilibrium, time-relaxation towards equilibrium) 
generically are exponentials; and their thermostatistical foundation is Boltzmann- Gibbs 
statistics (i.e., q = 1). But weakly ergo die /mixing phenomena also are ubiquitous in Nature 
(e.g., biological, socio-economical, human cognitive phenomena, etc); essentially, they are 
driven by microscopic interactions which are long-ranged in space time (long-range forces, 
nonmarkovian memory, fractal boundary conditions); their basic geometry tends to be dis- 
crete, multifractal-like; their thermodynamics is nonextensive; their central laws (energy 
distribution at equilibrium, time-relaxation towards equilibrium) generically are power-laws; 
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and the thermostatistical foundation of (at least some of) them (hopefully !) is the q ^ 1 
statistics. The allowance for nonextensivity, in general, and for a nonextensive entropy, in 
particular, appears to be the "price" to be payed in order that Boltzmann's "mechanical" 
(i.e., one system evolving along time) manner of thinking about macroscopic systems coin- 
cides, at the level of the concrete mathematical results to be compared with the experimental 
data, with Gibbs' "ensemble" (i.e., many systems at a fixed time) manner of thinking. This 
coincidence of results is, since one century, well known and understood for standard sytems. 
Our aim here is to extend it to a large variety of anomalous systems. 

The second of these conjectures is, at the present moment, so hard to rationalize that I 
dare to mention it here only because, after having been exposed to so many mathematical 
and physical arguments (that have been included in the present review), the reader might 
accept to honor me with his (her) indulgence, and have a look at the following few, intuitive 
lines. I believe that a deep analogy (maybe a kind of isomorphism, through the use of math- 
ematical structures like the co-homology groups) exists between crystallographic structures 
such as crystal - quasicrystal - fluid, and nonlinear dynamics such as integrable - (weak) 
chaotic - (strong) chaotic. In some sense, they appear as space and time versions of the 
same mathematical structures. The first case concerns crystals (i.e., <i-dimensional Bravais 
lattices) and integrable dynamics (i.e., motion on simple orbits), and its essential invariance 
is the discrete translational one. The third case concerns strongly disordered systems like 
fluids (liquid, gases) and strongly chaotic dynamics, and its essential invariance is the contin- 
uous translational one. Finally, the second, and intermediate, case is by far the most subtle 
one (and probably this is why it is the one that humanity took the longest time to discover), 
and concerns quasicrystals (e.g., Penrose tilings, amorphous substances like glasses, spin- 
glasses, and other structures known to have (multi)fractal scalings; probably most of the so 
called complex spatial phenomena belong to this group) and weakly ergodic dynamics (e.g., 
edge of chaos, strange attractors, self-organized criticality, probably most of the so called 
complex time phenomena); its essential invariance is the dilatation one. In the first case 
we have the (space or time) highest predictability, and statistical methods are out of place. 
In the third case we have the (space or time) lowest predictability, and statistical methods 
exhibit their full power. Finally, in the second case, we have an intermediate predictability, 
and the statistical methods have to be "intrinsically nonlinear" in some sense, in order to be 
applicable and useful. There will be no surprise for the reader if, at this point, I admit that 
I believe that the statistical mechanics to be associated with the third present case of course 
is the BG one, whereas it might be the q 7^ 1 statistical mechanics the one to be associated 
with the present second case ! (see also [|138|| , where some preliminary, but nevertheless 
concrete, calculations exhibit this kind of connections). Let us now remind that Wiles' 1995 
celebrated proof ||139|| of Fermat's last theorem was deeply related to quasicrystals since it 
was based on the proof of the Taniyama-Shimura conjecture about modular elliptic curves 
and used certain Hecke algebras. Consequently, I hardly dare to explicitly state a simple 
and unavoidable corollary, namely that, if my present second conjecture turns out to be, 
in some nontrivial and precise sense, correct, then the q 7^ 1 statistical mechanics must be 
related to Fermat's last theorem ! 

Through the complete analysis, in more detailed terms, of the various aspects tackled in 
the present review, we could learn a lot and, very especially, (precisely) when the celebrated 
Boltzmann factor is the correct theoretical description of natural systems at thermal equi- 
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librium ! This famous and so useful factor would then become, not a "dogma" , as referred 



to by Takens ||140|| , but a theorem ! 
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FIGURE CAPTIONS 

Figure 1 - W = 2 illustration of the escort probabilities: = pq+ ^l_ p y ■ 

Figure 2 - Generalization (Eq. (42)) of the Boltzmann factor (recovered for q = 1) 
as function of the energy E at a given renormalized temperature T', assuming a constant 
density of states. From top to bottom at low energies: q — 0, 1/4, 1/2, 2/3, 1, 3, oo (the 
vertical line at E/T' — 1 belongs to the limiting q = distribution; the q — > oo distribution 
collapses on the ordinate). All q > 1 curves have a (T" / ' E) q ^ q ^ lS) tail; all q < 1 curves have 
a cut-off at E/T' = 1/(1 - q). 

Figure 3 - Log-log plot of some cases like those of Fig. 2 (T" = 1, 5 for each value of q). 

Figure 4 - Central conjecture of the present work, assuming a Hamiltonian system 
which includes two-body (attractive) interactions which, at long distances, decay as r~ a . 
The crossover at t — r is expected to be slower than indicated in the figure (for space 
reasons) . 

Figure 5 - The one-jump distributions p q (x) for typical values of q. The q — > — oo 
distribution is the uniform one in the interval [— 1, 1]; q — 1 and q = 2 respectively correspond 
to Gaussian and Lorentzian distributions; the q — > 3 is completely flat. For q < 1 there is a 
cut-off at \x\/a = [(3 - q)/(l - q)} 1 ' 2 . 

Figure 6 - The g-dependence of the dimensionless diffusion coefficient A q (width of the 
properly scaled distribution p q (x,N) in the limit N — > oo). In the limits q — > 5/3 — and 
q -> 5/3 + we respectively have A q ~ [4/9]/ [(5/3) - q] and A q ~ [A/(9ir^ 2 ]/[q - (5/3)]; 
also, lim 9 ^ 3 A q = 2/ix 1 ! 2 . 

Figure 7 - The index q to be associated with a person whose critical value corresponding 
to Eq. (69) is V c . People with q < 1 (q > 1) tend to avoid (seek) risks for that particular 
game. The case q — 1 corresponds to an ideally rational agent. 

Figure 8 - Distribution of the transverse momenta Pt obtained in electron-positron 
frontal collisions of energy W varying from 14 to 161 Gev. The dotted line corresponds to 
q — 1 (i.e., a Hagedorn type of fitting as given by Eq. (87)) for all values of W. The solid 
lines correspond to q ^ 1 fittings. 
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Figure 9 - The values of q and T used in the fittings of Fig. 8. When W approaches zero, 
q approaches unity, i.e., Hagedorn's theory; T is essentially insensitive to W, as physically 
desirable. 

Figure 10 - Time evolution of £ = N(t)/N(0) associated with MbCO in glycerol-water. 
Dots: experimental data. Dashed lines: fittings with Frauenf elder's empiric law (Eq. (88) 
or Eq. (90)). Solid lines: fittings with the solutions of Eq. (91) (see Fig. 11). 

Figure 11 - Temperature dependences of the parameters used to fit the experimental 
data of Fig. 10. 

Figure 12 - Distribution of the "vertical" velocities during diffusion of endodermal 
Hydra cells in an ectodermal aggregate. The abcissa units are 10~ 6 m/hour. The fitting 
was obtained using q = 1.53 (see the text). 

Figure 13 - Distribution of ISI and PRD papers having received x citations, (a) and 
(6) exhibit the fittings in ||109|| ; (c) and (d) exhibit our present fittings (see the text). 

Figure 14 - Distributions of price changes for US dollar-German mark exchange rates 
and fittings using assymetric g-distributions (see the text). 

Figure 15 - Distributions of velocity differences and fitting using assymetric q- 
distributions (see the text). 
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TABLES 

TABLE I. Some useful formulae written, through Eqs. (9) and (10), in a Boltzmann-Gibbs-like 
form 



Equiprobability entropy 



S q = kln q W 



Generic entropy 



S q = -fc(lnqp)^ 



Canonical equilibrium distribution 



-f3(H-U q )/Trp-i 



Pq 



Tr e, 



0(n-U„)ITr 



re -d'n ~ Trp q + (l- 



q)f3U 



Z q = Tr e -M-U q )/Tr Pq { ^ ^ _ ^ 



Partition functions 



Internal energy 



U„ 



° In 7 

- m g Z/ g 



Free energy 



Fg — Ug ~ T.S'g — — g ln g 



Anomalous diffusion probability distribution 



p 9 (z) 



-/3a^ 



Agg 
Ax(0)^0 Ax(0 



Sensitivity to the initial conditions (d = 1) 



lim 



W q ({ Pl }) « e^ (tol) 



Likelihood function 



U(N,T) 



Power-law interactions (oc R a ) 



NN 



u(±) (N = N*[l + (a/d)]; N* ee ln [a/d] N) 



T(t) _ In, [TJW 



Simulated annealing (cooling rythm) 
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12 3 4 



E/r 



p(t,N,E) 



EXTENSIVE SYSTEMS ( a > d ) 



V 



- lim lim = lim lim 

N— oo t^oo t — oo N^oo 

(q = 1 ) 



T 



oc g(E) exp( -(3 E ) 
(N* not neede d) 



t 



NONEXTENSIVE SYSTEMS ( 0< a < d ) 



p(t,N,E) 



(q*1 )v 



— lim lim 

t — *oo N— »oo 



ocg(E) [1- (1-q) P'E] 



q/(1-q) 



t 



V 



— lim lim 

N— oo t^oo 

(q = 1 ) 



oc g(E) exp( -p E ) 
V [N* needed) 



X(N) t 



with lim T(N) =oo and 

N— oo 



N* = 



N 1-a/d -1 
1 - a/6 



2 -r 

q 

1 - 



o - 
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